Spatial and temporal turbulent velocity and vorticity power spectra from sound 

scattering. 
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By performing sound scattering measurements with a detector array consisting of 62 elements in 
a flow between two counter-rotating disks we obtain the energy and vorticity power spectra directly 
in both spatial and temporal domains. Fast accumulated statistics and large signal-to-noise ratio 
allow to get high quality data rather efi'ectively and to test scaling laws in details. 



PACS numbers: 43.30-Fm,43.35-Fd,47.32-y 

One of the challenging experimental tasks in studies of 
turbulent flows is developing new tools to measure spec- 
tral characteristics of velocity and vorticity fields in a spa- 
tial domain. Such measurements will allow direct com- 
parison of experimental data with a theory without ex- 
ploiting the Taylor hypothesis particularly in those cases 
when its use is rather questionable. Currently, from con- 
ventional methods in use only particle image velocimetry 
(PIV) is capable of performing this task on the expense of 
a very heavy data analysis to get representative statistics 
and a reasonable signal-to-noise ratio. 

In this Letter we present results on direct measure- 
ments of the velocity and vorticity fluctuations spatial 
and temporal spectra by using sound scattering ampli- 
tude measurements. We achieve this goal by simultane- 
ous acquisition of sound pulses on 62 sound detectors, 
arranged in front of a linear emitter in the same plane, 
as has been already described in Ref. [5] . There were sev- 
eral attempts in the past to experimentally probe turbu- 
lent flow characteristics by the sound scattering method. 
These experiments were mostly concentrated on studies 
of vorticity distribution [1] as well as temporal dynamics 
of vorticity fluctuations [2,3] and a large scale circulation 
in high Reynolds numbers flows [4,5]. 

As we pointed out [5] the main obstacles to get reli- 
able and complete information about the velocity field 
from sound scattering are a finite width of a sound beam 
and far field approximation required by the theory [6-8] . 
Since an integral scale of the velocity fluctuations is of 
the order of the beam width, the former limitation is re- 
leased for the turbulent velocity fluctuations. Thus, to 
get spatial information about the velocity spectra with 
sufficient dynamical range from the sound scattering data 
one needs to exploit a large number of elements in the 
sound detector array, to use lock-in detection technique 
to improve signal-to-noise ratio, and to apply the Huy- 
gens far-field construction method [5]. Then the ampli- 
tude of the complex wave function of the scattered sound 
in the far-field limit can be related to the Fourier trans- 
form of the velocity field (and also the vorticity) in a 



scattering plane as follows [7,^ 
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where ^scat = 4* — ^'rest, ^ and '^rest are the com- 
plex wave functions that describe the sound pressure os- 
cillations in the presence of a flow and without it, re- 
spectively; fco and c are the sound wave number and the 
velocity, respectively; 9 and kg are the scattering angle 
and the wave number, respectively, that are related to 
each other via kg = |fcs| = 2kos'm9/2; and r = |r| is 
the distance from the center of a scattering region till 
a detector. F^r {vx} is the 2D Fourier transform of the 
velocity component in the forward direction of the beam 
that is related to the Fourier transform of the vorticity 
via F,r^ {vx} = ^ cot (0/2)F,-; {(V xv),}. 

Since in many cases there exist severe technical obsta- 
cles to obtain a signal in a far-field limit, a construction of 
the far-field scattering wave function from the near-field 
measurements was suggested and used [5] . The method is 
based on a mathematical description of the Huygens prin- 
ciple in optics and derived from the Rayleigh-Sommerfeld 
integral [9] as: 
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where rd and rf are the distances measured from the 
cell center till the detector and the far-field region, re- 
spectively; 'i{0,y')scat and \E'(a:, are the scattering 
wave functions at the detector (as measured) and at the 
far-field, respectively [5]. 

The experimental set-up is described in details in Ref. 
[5]. The von Karman swirling water flow was pro- 
duced between two counter-rotating disks of a diameter 
2R = 280 mm with four triangular blades of 20 mm high 
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and 5 mm thick and with rims. They are driven by two 
dc brushless motors, which velocity is controlled with a 
stability of about 0.1% via optical encoders. This set-up 
is well recognized to generate a strong intensity turbulent 
flow in a confined region (see for example [2]). The flow is 
confined by a perspex cylinder of an inner diameter 290 
mm and 320 mm in height and disks separation 205 mm. 
By changing a rotation frequency the Reynolds number, 
Re = 2Q.R'^/v, is varied between 2.5 • 10^ and 1.7 • 10^ 
that corresponds to the Taylor microscale Reynolds num- 
ber Rx between 200 and 570. Here is the angular 
velocity of the disk, u is the kinematic viscosity, and 
the energy dissipation, e = 4.9 • 10^^*i?,e"^ W/kg, and 
the rms of the velocity fluctuations in the middle plane, 
Vrms = 0.5 • 10~^i?e m/sec, obtained from the global 
torque and the hot wire anemometry (HWA) measure- 
ments, respectively [10]. 

The sound scattering measurements were conducted in 
the middle plane between the disks and in the plane at 30 
mm below it. The emitter was 100 mm long and 10 mm 
wide, and the size of the scattering region was defined by 
the width of the detector designed as a linear array of 62 
acoustic detectors with 1 mm spacing and 62 x 10 mm 
active area (from Blatek). A range of frequencies cov- 
ered in the experiment was between 1 and 7 Mhz. The 
acquisition system is built in a heterodyne scheme that is 
based on 62 lock-in amplifiers with 62 preamplifiers. The 
details of the design and its operation are presented in 
Ref. [5]. A typical sound propagation time through the 
cell is about 200 /isec that is a typical freezing time seg- 
ment. Within this period one pulse is sent, and the flow 
is almost frozen. Since an each sound pulse has a suffi- 
ciently high signal-to-noise ratio, it is used to construct a 
sound far-field wave function from sound scattering sig- 
nals acquired by all detectors. Obtained as a function 
of time the wave function provides a possibility to study 
statistics of velocity (vorticity) fluctuations in a turbu- 
lent flow. 

A scattering wave function emitted from a single trans- 
ducer and obtained by the detector array provides infor- 
mation about the velocity (vorticity) field only on a sin- 
gle curve in a {{ka)xT{ks)y) plane (sec Fig.l, left inset, 
solid line). To get information on other curves in the 
plane (see Fig.l, left inset, dashed fines) requires to use 
sound beams emitted by different transducers in many 
different directions simultaneously. Then complete struc- 
ture functions of the velocity (vorticity) fluctuations can 
be retrieved without any assumption about isotropy and 
homogeneity of a turbulent flow. However, with only 
one emitter and the detector array we should rely on 
the isotropy and homogeneity assumptions of the tur- 
bulent flow under studies. In this case the sound scat- 
tering provides direct measurements of the energy spec- 
trum as well as the Fourier transform of the vorticity 
structure function in a spatial domain that are related 



as E{ks) = ^-\F^ {(V X 'y)z}p. Here the kinetic energy 
per unit mass is deflned as / E{k)dk = ^ J f^'^^?', where 
A is the cross-section area of a sound beam. We neglect 
velocity variations in the direction perpendicular to the 
beam propagation by averaging over the beam thickness. 
In Fig.l we present a typical result on the time-averaged 
velocity Fourier transform in the far-field obtained from 
the scattering wave function via Eq.(l) and the far-fleld 
reconstruction according to Eq.(2). The same function 
but observed directly in the detector plane looks drasti- 
cally different (see right inset in Fig.l). These data are 
the result of averaging on 60,000 sound pulses at 1.8 kHz 
repetition rate and at frequency 3MHz taken in the von 
Karman swirling flow at Re = 1.5 • 10^. 

In Fig. 2 the resulting energy spectrum as a function of 
the wave number is shown. The dotted line denotes the 
"-5/3" slope according to the Kolmogorov law to demon- 
strate that indeed in some range of the wave numbers 
the spectrum follows the expected law. It can be com- 
pared with the results on the energy spectra obtained by 
PIV, where though a shorter scaling region is also ob- 
served [10]. In the sound measurements shown in Fig.2 
the wave numbers are limited by the range of values 
0.1 < k < 0.8 mm~^, outside of which the spectrum 
cannot be retrieved. The lower side of the range of k is 
limited by the beam width, i.e. one cannot get informa- 
tion on a scale exceeding the beam width (in the energy 
spectra from PIV the scaling region begins already at 
k = 0.06 mm~^ for Re > 10^) . On the higher side of 
the range the ultimate limit is determined by the size 
of the element in the detector array, namely kmax < 27r 
mm""'^ corresponding to 1 mm detector array spacing. 
However, even before this limit is reached, the energy 
spectrum is cut on the higher wave number side by the 
visibility at large scattering angles through the detector 
aperture. It means that some sound rays are blocked by 
the limited length of the detector array. According to 
the detector array length and the cell diameter at an- 
gles exceeding about 6° the visibility starts to deterio- 
rate (see the inset in Fig.2). The corresponding limit is 
^S(Hm) ^ 2fco sin ((?ma2:/2) — O.lfco- It gives about 1.3 
mm~^ at 3MHz compared to 0.8 m^^ observed. We 
tested this relation at various frequencies, and the results 
on the energy spectra taken at different sound frequen- 
cies are presented in Fig. 3. One finds that by changing 
the frequency from 1 up to 7MHz the upper wave number 
limit moves linearly toward the highest value of about 1.5 
mm~^, but still far away from kmax = 27r mm"-"^. 

Useful information in terms of the turbulent flow en- 
ergy dissipation, e, can be gained from combined pre- 
sentation of the energy spectra obtained from the sound 
scattering measurements at different Reynolds numbers. 
There is a known scaling law for the energy spec- 
tra [11-13] that appears as the result of plotting the 
scaled energy density spectrum e~'^/^r]~^/^E{k) versus 
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kr], where 77 is the Kolmogorov dissipation scale defined as 
r] = (i^^/e)^/^ [13]. The idea is to find the best match be- 
tween the scaled spectra at different Re with fitting val- 
ues of e. It turns out that the scaling exists for the data at 
all values of Re, and the results for the sound scattering 
in the middle plane (curve 1) and in the plane at 30 mm 
below it (curve 2) arc shown in Fig. 4. Each data set for 
each Re consists of 4 • 10® points. The dependence of the 
energy dissipation, e, on Re is found to be e ~ i?e^'^^'^'^ 
with a good quality of the fits. The exponent value is 
rather close to 3, the expected one according to the di- 
mensional analysis [10,14] (compare with the expression 
for e presented above). The Kolmogorov constant C in 
the Kolmogorov equation e~'^^^ri~^^^E{k) = C{kvi)~^/^ 
is determined experimentally from the fit of the plots in 
Fig. 4, and the value is C ~ 0.8 for the curve (1) and 
C ~ 0.9 for the curve (2). 

We also calculate the integral scale in the 
flow at Re = 1.2 • 10*^, using [15] L,„t = 
(37r/4) / k-^E{k)dk/ J E{k)dk = 40 ± 10 mm, based 
on the system scale 0.02 mm~^ and the scaling region 
0.06 < k < 1.5 mm""'^. This value occurs to be rather 
close to the beam width, and the corresponding wave 
number is located close to the lower end of the wave 
number range of the spectrum. 

To test the Taylor hypothesis [13,15] for the swirling 
flow the energy spectrum in the frequency domain is cal- 
culated as E{f) = i E{t) e^{-i2'Kft)dt and E{t) = 

I ■IF'l-^fe*'-^ ^ v{t))z\'^d'^ks. A proper energy spectrum 
can be obtained only, if the lowest ks are available 
due to the pole at k~^ in the integrand. So we turn 
out to the Fourier transform of the vertical vorticity. 
In Fig. 5 (upper curve) we present the time-averaged 
power spectrum of the vorticity as a function of k 
(based on 2 • 10^ points) at Re = 1.2 • 10® obtained at 
the sound frequencies 2.5 MHz (asterisks) and 5.8MHz 
(diamonds). The dashed line with the scaling expo- 
nent "-2/3" represents the expected dependence accord- 
ing to the Kolmogorov predictions [13,15]. It can be 
compared with the space-averaged over the beam area 
Fourier transform of the enstrophy in the frequency do- 
main N{f) = i N{t)eyip{-i2Trft)dt, where N{t) = 

\U{Vxv(t)U^d^r = {2n)'J\F^-^{Vxvm^d^ks. 
The latter shows a rather wide scaling region (see the 
lower curve in Fig. 5). The advection velocity, Vt, found 
from the two plots in Fig. 5 is Vr = 2nf /k ~ 0.6 m/sec 
that is rather close to the average velocity in this exper- 
iment [10]. 

We would like to point out also that the mean en- 
strophy value, marked on Fig. 5 by asterisk on the left 
hand-side ordinate axis, is about an order of magnitude 
larger than one corresponding to a rigid body rotation 
with the same rotation velocity of 300 rpm and the vor- 
ticity of 63 sec~^, or the enstrohy of about 4000 sec~^. In 
the co-rotational disks geometry at the same parameters 



and with a large scale single vortex flow configuration the 
enstrophy is also much lower [5]. 
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FIG. 1. Time averaged velocity structure function obtained 
by sound scattering via the far-field construction taken at 

Re = 1.5 ■ 10® and 3 MHz. The right inset: the same property 
obtained from sound scattering data at the detector without 
far-field construction. The left inset: path in the wave num- 
ber piano (solid line), on which information on the velocity 
(vorticity) spectrum for a given beam direction is obtained; 
other curves (dashed lines) contain information on scattering 
from sound beams emitted in various directions. 



FIG. 2. Energy spectrum derived from the data presented 
in Fig.l. The inset: scheme of the aperture limit for the sound 
detector array. 



FIG. 3. (color online) Energy spectra taken at various 
sound frequencies from 1 up to 7 MHz at iie = 1.2 • 10®. 
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FIG. 4. (color online) Scaled energy spectra as a function 
of reduced wave number at various Re: (1) at h = (2) at 
h = —30 mm. The dashed lines show the Kolmogorov scaling 
law with the exponent "-5/3". 



FIG. 5. Upper curve: Power spectra of time averaged vor- 
ticity as a function of k taken at J?e = 1.2 • 10^ and frequencies 

2.5 (asterisks) and 5.8 (diamonds) MHz. Lower curve: Power 
spectra of space-averaged enstrophy as a function of /. The 
dashed lines show the Kolmogorov scaling with the exponent 
"-2/3". 
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